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.  ,/J  %  INTRODUCTION 

Since  the  linear  solution^  for  the  problems  of  f'  ee  surface 
waves  were  found  by  Kelvin  (T.887)  and  Michell  (1898),  theoretical 
and  experimental;;  regeargh  on  the  improvement  of  ship  forms  by  the 
application  c.f  the  linear  theory  has  been  conducted  by  many  naval 
hydrodynam.ic.1sts,  and  many  fruitful  results  have  been  obtained. 
However^  the;  linear  theory,  particularly  for  surface  ship  problems, 
seems  ta.  ha  ye  severe  limitations.  Yet,  the  higher  order  wave 
theory  of  surface  ships',  being  inherently  difficult,  has  developed 
at  an  ex t>  eme& y  0 w  pace.  Recently,  because  of  the  easier  avail¬ 

ability  and  improvement  of  the  computing  speed  of  high  speed  com¬ 
puting  machines,  this  complicated  by  highly  necessary  theory  has 
become  more  and  more  attractive.  Wehausen  (1963)  considered  the 
systematic  development  of  higher  order  ship  wave  theory.  Sesov 
(1961)  formulated  the  second  order  wave  resistance.  Yim  (1966), 
considering  the  free  surface  condition  more  carefully,  investi¬ 
gated  the  higher  order  potential  and  the  wave  resistance.  This 
author  also  investigated  the  effect  of  the  line  Integral  on  the 
free  surface,  among  many  other  higher  order  effects  (Yim,  1964), 
unfortunately  with  some  computational  error. 

In  this  report  another  attempt  on  the  investigation  of  the 
line  integral  on  the  free  surface  is  now  made,  with  8  parabolic 
sh.'^  form  and  by  a  much  simpler  method  of  computation.  As  before, 
the  influence  of  the  line  integral  on  the  asymptotic  wave  height 


is  calculated. 
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It  has  lone  been  known  that  a  pressure  point  moving  on 
the  free  surface  produces  a  peculiar  flow  field  along  its  path, 
Lamb  (19^5),  Ursell  (i960}.  Although  a  pressure  point  is  known 
to  be  equivalent  to  a  doublet  point,  not  much  information  is 
available  about  the  properties  of  general  singularities  on  the 
free  surface.  However,  due  to  the  recent  development  of  slender 
ship  theory,  Vossers  (1962),  Tuck  (1963)*  Maruo  (3. 962),  Joosen 
(1963),  and  the  higher  order  ship  wave  theory,  the  importance 
of  the  singularity  on  the  free  surface  increased  rapidly.  Thus 
the  author  has  devoted  his  attention  to  obtain  as  much  infor¬ 
mation  as  possible  abouc  the  behavior  of  singularities  on  the 
free  surface  -  before  calculating  the  line  integral  due  to  a 
singularity  distribution  along  a  line  on  the  free  surface. 

SINGULARITIES  ON  THE  FREE  SURFACE 

We  locate  an  origin  of  the  rectangular  right-handed  co¬ 
ordinate  system  on  the  undisturbed  free  surface  (or  a  mean  free 
surface)  with  the  directions  of  three  axes  as  shown  in  Fig¬ 
ure  l.  If  we  consider  a  pressure  point  at  the  origin  moving 
with  constant  velocity  -V  on  the  free  surface,  the  flow  field 
satisfies  Laplaces  equation 

7s*  «  0  [1] 

with  the  free  surface  boundary  condition 

*xx  +  Vx  +  *  0 

-  k/v*)  1?) 
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on  z  =  0,  except  at  the  singularity,  and  with  the  proper 
boundary  conditions  a t  infinity.  The  solution  is  well  known 
(Wehausen,  I960;  Peters, 1948;  Ursell, i960). 
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[3] 

This  Is  exactly  the  same  as  the  limiting  case  -'f  zero  submer¬ 
gence  of  the  potential  due  to  a  point  doublet  with  its  direc¬ 
tion  -V. 

The  pressure  point  being  naturally  the  limiting  case  of 
zero  submergence  of  the  lift  element,  this  fact  can  be  seen 
easily  from  the  free  surfac  boundary  condition  outside  the 
singularity, 

cp  +  k  ®  »  0  on  z  ■  0 

XX  o  z 
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If  we  consider  a  point  singularity  at  a  point  (xi ,  yi ,  0)  and 
the  flow  which  satisfies  [4]  at  a  point  (x,  y,  0)  we  can  re¬ 
write  [4] 

«P  v  (xx  *yi  ,0;x,y,0)  +  k  ®  (xx  ,yx  ,0;x,y,0)  =  0  [5] 

XI  Xi  O  Zi 

since 

m  =  -cp  ,  <P„  =  -®  ,  cp,  *  +<P„  C6] 

Xx  x  yx  y  zx  2 

referring  to  the  Green's  function  (see, e .g., Lunde  (1951))  which 
represents  a  potential  due  to  a  point  source.  Then  from  [5! 

• 

<PV  (0,yx  ,0jx,y,0)  -  k  f  ( p  (xx  ,y4 ,0;x,y,0)dXx  [73 

Xx  o  J  Zx 

0 

This  means  that  the  flow  on  z  ■  0  due  to  a  point  doublet  located 
on  the  free  surface  is  the  same  everywhere  on  z  ■  0  as  that  due 

to  a  vortex  element  on  the  free  surface.  Thus  the  two  flows  are 

identical  everywhere  in  the  fluid. 

A  point  source  located  on  the  free  surface  can  also  be  con¬ 
sidered  in  a  similar  manner  as  the  point  doublet  on  the  free 
surface.  Since  a  j>int  source  at  the  origin  may  be  considered 
as  a  line  distrl  ion  of  doublets  on  the  positive  x  axis,  a 
point  source  on  the  free  Surface  can  be  interpreted  as  the  uni¬ 
form  pressure  distribution  along  the  x  axis. 
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The  wave  height  C  due  to  the  pressure  point  is,  according 
to  the  dynamic  relation, 

•  •  > 

o 

and  Equation  [3]  with  some  changes  of  variables,  can  be  ex¬ 
pressed  (see  Ursell  (i960))  as, 

C(x,y)  »  Ci  (x,y)  +  Ca  (x,y) 


•  ■,|„i  _  1/2 

-1—  fd K  K»e  f dm  3ln  ro  coafKyc.osm) 

x^pV*  J  J  kQa  +  It  sin*m 

o  o 


x  exp  {.-k0*(i  >  u*  ))  du  C  93 

where  Ci  la  the  local  disturbance  which  decays  rapidly  with 
increasing  | x| ,  and  Ca  la  the  regular wave*  The  regular  wave 
Ca  far  behind  the  pressure  point  can  be  decomposed  into  two 
par‘  m  the  transverse  wave  and  the  divergent  wave. 
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If  the  pressure  is  distributed  on  a  finite  area,  there 
exists  no  singularity  in  the  velocity  field  (Stoker,  195-7,  or 
Ursell,  i960).  However,  Equation  C 93  has  an  innate  singular¬ 
ity  at  the  pressure  point  and  near  its  path  such  that  the  wave¬ 
length  is  Infinitely  decreasing  and  the  wave  height  is  in¬ 
finitely  increasing  as  y  -  0  (Lamb, 19^5,  Ursell, i960),  although 
the  wave  height  on  y  =  0  or  the  flow  in  2  <  0  is  finite.  The 
singularity  near  the  path  arises  from  only  the  divergent  wave, 
while  the  transverse  wave  has  a  finite  amplitude  except  at  the 
origin.  Tne  wave  resistance  due  to  the  pressure  point  also 
blows  up.  If  we  follow  exactly  the  same  procedure  as  in 
Ursell* a  paper  (i960),  it  is  obvious  that  the  flew  field  due 
to  a  point  source  on  the  free  surface  has  exactly  the  same  kind 
of  singular  phenomena  with  a  slightly  weaker  singularity.  In 
general,  a  line  doublet  distribution  on  the  free  surface  also 
has  the  same  kind  of  singular  phenomena ,  If  the  distribution 
of  doublet  strength  near  the  end  points  smooths  sufficiently 
to  Eero,  the  flow  field  be  haver  nicely,  distrlbu- 

.  %  lor- ;  along  a  straight  line  parallelte  %  In  general  has  neither 
finite  wave  resistance  nor  finite  wave  height  oh  the  path,  un¬ 
less  the  end  distribution  is  smooth  enough.  However,  if  the 
line  wnere  the  sources  are  distributed  is  smooth  and  not  paral¬ 
lel  to  V,  the  flow  may  behave  well.  These  results  mentioned 
above  can  be  easily  obtained,  if  we  apply  the  method  of  steepest 
descents  to  each  expression  using  the  result  in  the  case  of  the 
pressure  point  In  Ureeil's  paper  (I960)  and  having  Havelock's 
formula  for  the  wave  resistance  (Have  lock, 193^)  In  mind. 
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The  singular  behavior  in  the  actual  flow  never  arises. 
Even  if  we  just  consider  the  surface  tension  (Wehausen,1963), 
the  divergent  wave  heights  have  no  singular  behavior  anywhere. 
Thus,  in  the  case  of  point  singularities  on  the  free  surface 
considered  above,  the  wave  height  near  the  path  of  singular¬ 
ities  may  well  be  considered  as  smooth  to  the  finite  theoreti¬ 
cal  value  at  y  =  0. 

FLOW  FIELD  DUE  TO  A  SHIP  AND  THE  LINE  INTEGRAL 

From  the  Green's  formula,  the  potential  cp  can  be  written 
using  the  Green's  function  G(§  ,T),C  ;x,y,  z )  which  satisfies 


on  C  =  0 


Outside  the  ship, 


in  z  ^  0 


where 

r8  *  (x  -  ?  )8  +  (y  -  h)8  +(z  -  C  i* 


[10] 


Cll] 


G  -*  0  for  |  r|  -•  *  in  z  <  0 


and 
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cp  =  -L  JJ C‘P(5,r»,C)Gn(5,ri,C,x,yJz)  -  dS  [  12] 
S 

where  n  is  the  normal  into  the  fluid,  and  S  consists  of  the 
ship  surface  Sg  and  the  free  surface  since  the  integral 
from  infinity  is  null. 

From  the  free  surface  integral. 


where  l  is  the  intersection  of  the  ship  surface  and  the  z  =  0 
plane.  This  line  integral  on  the  free  surface  has  been 
unusually  neglected  in  Michell's  ship  theory.  A  numerical 
computation  of  the  effect  from  this  line  integral  to  the  bow 
wave  on  the  centerplane  far  behind  the  wedge  ship  was  con¬ 
sidered  (Yim, 1964).  Now  we  consider  a  parabolic  forebody  with 
a  parallel  midsection  represented  by  a  source  distribution  on 
y  ■  0,  z  s  o 


■i 
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Slnce 


f=ao?  in  0  <  §  <  1 


G~  =  -  G 


x 


we  have 


T  - 


“x  2irk 


w*  Sr]/- 2  f  Vs  Sr + 

o  o 


f  daf 

J  §  d?a 


d§ 


27Tk 


-  <p(.0)Gx(x,0)  -  2  /  Gx«P?(l  d? 


-  G(x,l)PP(l)  +  G(x,0)cp(0)  +  /  GCp?d? 


C 183 


neglecting  the  higher  order  terms.  The  first  term  represents 
the  wave  height  due  to  a  point  source.  This  is  singular  as  we 
have  seen  In  our  previous  section.  This  phenomenon  is  familiar 
to  hydrodynamic ists  especially  concerned  with  higher  order  prob 
lems,  The  singularity  usually  becomes  more  severe  when  the 
order  becomes  higher.  Thus,  the  reliability  of  the  higher 
ordex  theory  in  most  problems  in  hydrodynamics  may  be  judged 
through  experiments  only. 
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We  consider  here  the  wave  height  on  y  -  0  far  behind  the 
ship  where  the  theoretical  wave  height  has  a  definite  value  if 
we  consider  that  the  divergent  waves  are  continuous  everywhere 
because  of  surface  tension  or  for  other  reasons. 

Using  Wigley's  notations  (1931)  and  Havelock’s  example 
(1932),  we  can  write  the  total  wave  height  in  0  <  x  <  1 


C  -  7— -  Q  (k  x)  +  Q  (k  ai  -  x  +  l) 
k  l  o v  o  ‘  o x  0 

O  \ 


-  jp  foi  (K  x)"  +  Qx  (k  1  -  x)  -  Qt  (k  ax  -  x) 

0  L 

+  Sfc.  (koiT^Ti  )]  -  4  [p0(kox)  -  i  p;1  <v)j 


where 


Qo(a)  *  I  lHo(t)  ’  Yo(t}}  dt 


H  (t)  :  Struve  function  of  order  zero 
o 

Y  (t)  :  Neumann  function  of  order  zero 
o x  ' 


r  a 

Qi  (a)  *J  Q.o(u)  du 

0  ra 

Po(a)  ■  -  I  J  yo(t) 


[19] 


’o 


P0(t)  dt 


C  203 
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This  value  of  C  is  plotted  in  Figure  2,  In  x  <  0 


2a 


C  =  - 


o 


Vkox'>  +  ir(Q‘<V>  -  +  nj 


+  Q  (k  ai 
ov  o 


x  +  1)  +  (kQai 


-  x)  -  Qi  (kQai  -x+l)j 


[21] 


where  x '  -  -x 

The  potential  is  unique  with  the  condition  at  infinity  cp  -*  0 
which  is  necessary  to  get  the  Green's  formula  (see,e .g. ,Kellog, 
1953 ) *  and  at  the  origin 


-«  -00 


Qi  (kQax 


+  1 )  + 


t/ 


Qi  (u)du 


l  f  0 

+  /  Qi  (u)du 

*o~l  - 


[22] 


Similarly 


on  z  ■  0 


[23] 


can  be  obtained. 
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*  iJ[po(kox)  -  ir{v'(v>  -  V'V  -  x>} 

-  P  [w>  *  i r  p  <V>  -  *  (ko~>}]  t26] 


By  application  of  the  method  of  stationary  phase  and 
neglecting  the  local  disturbance 


«  A;  sin  |kQx  +  jp  +  8  j 

for  x  >  >  1  [273 
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Then  the  total  wave  height  Ct  will  be 

C*  *  Cl  +  c£ 

=  Ai  sln|kQx  +  ^  +  S  J  +  Aa  sin|kQx  +  ^  +  /3 
=  AaSln|kox  +  jp  +  8  +  &J 

where 


A»  *  {Ai*  +  Aa*  +  2  A^Aa  cos  (0  -  6)}^ 

As  and  6  are  plotted  in  Figures  3-5* 

DISCUSSION 

The  calculated  result  shows  that  the  line  integral  effect 
on  the  wave  heignt  due  to  the  fore  body  of  the  parabolic  ship 
with  parallel  mid  body  appears  as  an  increment  of  the  ampli¬ 
tude  ard  a  forward  shift  of  the  wave  phase.  The  oercentage  Is 
more  in  lower  Froude  numbers.  Thus,  the  discrepancy  between 
the  experiment  and  theory  may  be  ascribed  to  higher  order  ef¬ 
fects  and  especially  to  this  line  integral .-effect .  Tft*  impor¬ 
tance  of  this  effect  can  also  be  shown  by  the  following  approach. 
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We  consider  a  double  model  ship  with  uniform  speed  -V  In 
the  infinite  medium,  where  z  <  0  correspond  to  the  case  of  zero 
Froude  number.  Then  the  pressure  Pi  on  z  =  0  will  be 

Pi  58  pf(<Pi  )  ~  Vpcpix 

where  p  is  the  density  of  water  and  cpi  is  the  potential  due  to 
the  double  model  ship.  Here  we  consider  on  z  *  0,  Kelvin's 
pressure  source  distribution  -Pi  which  satisfies 

^ x 

<Paxx  +  kQ9»z  -  -  -pf  on  z  «  0 

where  <p»  is  the  potential  In  z  *  0  due  to  Kelvin's  pressure 
distribution.  If  we  consider  the  ship  is  symmetric  with  re¬ 
spect  to  y  m  0  plane,  ft  la  also  symmetric  with  respect  to 
x  axis.  Thus 

9„  *  0  Oh  y  ■  0  V. 

w  .  ■ 

If  we  superpose  the  two  fldi^ 

9  *  9i  ♦  9t 

will  exactly  satisfy  Laplace's  Equation  Cl],  the  free  surface 
condition  on  z  «  0  and  the  linear  boundary  condition  of  the 
thin  ship  surface  on  y  -  0.  Thus  9  will  be  exactly  the  sane 
as  the  Mlcnell's  solution. 
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If  we  want  to  consider  the  problem  more  accurately,  we  have 
to  first  notice  the  fact  that  the  free  surface  is  limited  to  the 
region  outside  of  the  ship.  The  pressure  inside  the  double 
model  ship  on  z  »  0  is  most  strong,  and  the  effect  on  the  ship 
surface  due  to  the  Kelvin's  pressures  located  outside  of  the 
ship  is  not  too  large  because  of  the  symmetry  of  the  pressure 
distribution  with  respect  to  y  =  0.  Thus,  Michell's  theory 
which  neglects  the  strong  extra  pressure  inside  the  ship  can  be 
expected  to  lead  tu  a  sizable  error  on  that  account. 
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FIGURE  2-  FIRST  ORDER  WAVE  HEIGHT  OF  PARABOLIC  SHIP 
WITH  PARALLEL  MIDDLE  BODY 


PARABOLIC  SHIP.  V/\TflT  =0.258 ,  x/L=  15. 


FIGURE  4-  INFLUENCE  OF  HIGHER  ORDER  WAVE  ON  THE  FIRST  ORDER  WAVE  FAR  AFT  OF 

PARABOLIC  SHIP  2/ro.  =  0.2,  x/L  =  15. 
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INFLUENCE  OF  HIGHER  ORDER  WAVE  ON  THE  FIRST  ORDER  WAVE  FAR  AFT  OF 
PARABOLIC  SHIP.  V/VgT  =  0.408, x/L=  15. 
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